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Abstract The quantum-classical transitions of the escape rates in a 
uniaxial spin model relevant to the molecular magnet Mn^Ac and a biaxial 
anisotropic ferromagnetic particle are investigated by applying the periodic 
instanton method. The effective free energies are expanded around the top of 
the potential barrier in analogy to Landau theory of phase transitions. We 
show that the first-order transitions occur below the critical external magnetic 
field h x = \ for the uniaxial spin model and beyond the critical anisotropy 
constant ratio A = \ for the biaxial ferromagnetic grains, which are in good 
agreement with earlier works. 
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Macroscopic quantum tunneling(MQT) of magnetization was intensively studied both 
theoretically and experimentally in the last decade Q. In recent years much attention 
has been attracted to the phase transition problem in quantum tunneling of magnetization 
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]. The interest to this problem arises mainly from the successful experiments on the 
molecular magnet Mn 12 Ac |5|-|7j . It has been shown that the quantum-classical transition 
of the escape rate is analogous to the phase transition and the general conditions for the 
first- and second-order transitions are also analyzed. Quite recently an effective free energy 
F = acf) 2 + 60 4 + c0 6 for the transitions of a spin system was introduced as in Landau 
theory of phase transition. Here a = corresponds to the quantum-classical transition and 
b = to the boundary between first- and second-order transitions. In the common sense the 
first-order transitions are difficult to find in real systems and a uniaxial spin model is one 
of the very few examples which would exhibit the first-order transition. Various theoretical 
methods were used to deal with the spin tunneling problem @-|13j- We present in this 



letter the periodic instanton calculations for the quantum-classical transitions of the escape 
rates for the uniaxial spin model and a biaxial anisotropy ferromagnetic particle. The latter 
provides another candidate in which we may observe the first-order transition at the proper 
anisotropy constant. 

A rather simple and experimentally important system which may exhibit the first-order 
transition is the uniaxial spin model in a magnetic field parallel to x-axis H x described by 
the Hamiltonian 

n = -DSl - H X S X (1) 

which is generic for problems of spin tunneling. This model is believed to have relevance 
to the molecular magnet Mn 12 Ac with D the anisotropy constant. It was found recently 
that there exists a critical value of the external field below which the first-order escape 
rate transition occurs 0. Here the periodic instanton method is used to deal with this 
problem which should be more accurate as we need not resort to the double well potential 
approximation. Using the method of mapping the spin model onto a particle problem [|3],|9] , 
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the equivalent particle Hamiltonian is TC = ^— + U(x), where 

U(x) = S(S+ l)D(h x cosh x- if (2) 
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and m = 1/2D, h x = H X /2SD, S = S + 1/2. The minimum of the effective potential, 
xq = cosh -1 (1 / h x ) , has been moved to zero and S ^> 1 is used throughout. Recently the 
single kink and kink lattice solutions for a class of quasi-exactly solvable potential model 
in field theory, the Double sinh- Gordon (DSHG) potential ||14|| , have been found and the 
statistical mechanics of DSHG kinks are studied. Our motivation is based on that the 
potential (||]) is a special case of DSHG model and corresponding instanton and periodic 
instanton configurations may be found following the similar procedure. 

The vacuum instanton solution is nothing but the zero-energy solution of the equation 
of motion in imaginary time r 



1 ( dx 
—m — 

2 \ dr 



- U{x) = 0. 



(3) 



The instanton/anti- instanton located at tq is given by 
x c (t) = ±2 tanh -1 



tanh — tanh " 
2 



1 



where tanh ^ 



l-h x 
1+fta 



and the corresponding Euclidean action is 



1 - hlSD 



(4) 



Sr 



+OC / 1 + Jl - h 2 x 

mx c dr = 2S In — 

-DO \ tlx 



l-hl\. 



(5) 



To understand the finite temperature tunneling behavior we construct the periodic in- 
stanton configuration satisfying the periodic boundary condition similar to the kink lattice 
method in double sinh-Gordon(DSHG) theory. This solution is obtained by integrating the 
equation of motion 
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\ m \Tr) ~ U{X) 



■E 



(6) 



and is given by 



x p (r) = ±2 tanh 1 



tanhxi sn ( — — — -, k 



(7) 



with 
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^E' — h x 
1 



tanhxi 
tanhx2 



(8) 



k 



2h x SD sinh x\ cosh x 2 



SD^l-[h x -VW 

where sn (r, k) is the Jacobi elliptic function with modulus k and the complementary modulus 



k! = vl — A; 2 . Equivalent ly, 



1 - ^ T v^ 7 _ 1 - K (1 + k 2 )h x =F A;' 2 (l + ^) + y^fc 2 + #4 

tann xx? — t== — ; ; (9) 

1 + 1 + h x (i + p)^. zp jfe/2(x _ + v /4/ i ^ 2 + fc' 4 

and the characteristic length of the periodic instanton 

/ 
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„ ( v l + k 2 )h 2 x - k' 2 (l - h 2 x ) + h x ^h 2 x k 2 + k' 4 



2\ -1 



V 



(10) 



[l + k 2 )h 2 x + h x ^Ah 2 x k 2 + k' A 

This description corresponds to the movement of a pseudo-particle in the inverted potential 
—U (x) with energy —E and E' = E/S 2 D. The periodicity of the solution (g) is t p (E) = 4/3, 
f3 = K(k)£p, where K(k) is the complete elliptic integral of the first kind. The topological 
charge of periodic instanton 



Q p = 2x p (K(k)) = 2x 1 



(11) 



is smaller than the vacuum instanton case 

Q- 



00 dxJr) 



dr 



dx = 2xn. 
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Similarly, the periodic instanton size £p is also smaller than the zero-energy one £. The 
potential and the instanton configurations are depicted in Fig. 1. The particle starts from 
the turning point —x\ at imaginary time —f3 and reaches the other turning points x\ at 
(3. After the same interval, at time 2/3, the particle returns to its original position, i.e., it 
tunnels through the barrier twice in the whole period. 

The Euclidean action of the periodic instanton configuration is 



mx\ l+E dr = W + 2E(3 



(13) 



where 



W 



a 4 - k 2 ) U(a 2 , k) + k 2 K{k) + a 2 (K(k) - E(fc)) ] 



(14) 



D£ P a 2 1 

Here a 2 = tanh 2 a;i < k 2 and E(fc), IT(a 2 , fc) are the complete elliptic integral of second and 
third kind, respectively. 

The temperature dependent tunneling rate can be estimated by a Boltzmann average 
over the tunneling probabilities from excited states with energy E which is approximated by 
the semiclassical WKB exponents, r„ ~ e ~ 2W ( E )_ T n fa e quasi-classical approximation the 
transition rate is given by T ~ e - F mm/T^ w he r e _p min is the minimum of the effective "free 
energy" M 



F = E + TW'(E), W'(E) = 2W(E) 



(15) 



with respect to E. In the case of second-order transition the crossover temperature is given 
by 



T (2) _ ^0 



7T 



(16) 



where u = y\U"(0)\/m is the instanton frequency. It is convenient to introduce dimen- 
sionless temperature and energy variables 







T 



T, 



(2) 



P 



AU -E 

AU ' 
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where AU = U ma _ x — U min is the barrier height. To investigate the phase transition behavior, 
we need to expand the free energy around the top of the potential barrier. Near the potential 
maximum^ ~ 0) the expansion of elliptic integrals up to order of k 6 is seen to be flT5] 



K(Jfe) 



7T 



E(*) = \ 



1 + k l + l k 4 + —k e + • • 
4 64 256 

I — — — —k 4 - —k e - ■ ■ 

4 64 256 



(18) 



n^anh 2 ^,^ = ^ + l(3-2h x )k 2 + -^-(32hl-8h 2 x -G0h x + 45) k 4 



71 



512 



(-256/^ + QAh 4 + 480h 3 x - 88h 2 x - 350h x + 175) k 6 . 
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The other parameters in the free energy, such as a 2 and £ p , are also calculated in the same 
way and we obtain the result 

F/AU = l + 4h x (d- 1) k 2 + Ah x [h 2 x + 2h x -l-9 [h 2 x + h x ~l))k 4 

+4h x [e [2hi + 3h 3 x - hi - ^h x + g) - (2hi + 4/£ - 4h x + l)) k 6 . (19) 

There exists an exact relation between k 2 and P 



k2 = -2h x -p + h x p + 2h x yr^p 

-2h x - P + h x P - 2h xV / T^P' [ ' 

Expressing k 2 in power series of P, we have 

F/AU = 1 + (9 - 1)P + -(1 - ^-)P 2 + — (1 - ^- + ^tt7)P 3 + 0(P 4 ) (21) 
7 v ; 8 V Ah x ! 64 v 3/i x 16/^2 ; v ; y ! 

which coincides with Ref. [0 exactly. There indeed exists a phase boundary between the 
first- and second-order transitions, i.e., h x — \, at which the factor in front of P 2 changes 
the sign. We conclude that this critical boundary is inherent in the DSHG model and plays 
dominant role in the tunneling process of the uniaxial spin system with an external magnetic 
field, eq. (1). 

Turning now to the computation of level splittings of excited states. We have a more 
general formula for the double- well-like potentials in WKB approximation |l6|Jl7| j3|1 



AE = ^--Uxp[-iy] (22) 

where ui{E) = is the energy-dependent frequency and t(E) is the period of the real-time 
oscillation in the potential well. This level splitting formula is valid for entire energy region 
< E < AU. The calculation of the period t(E), equivalently the normalization constant 



of WKB wave functions [1181 , results in 



, dx 2K(k') 
t(E) = V2m / . = v ; . (23) 

Cl " W) SDjl-(VW-h T ) 2 



For energies near the bottom of the potential well the energy-dependent frequency tends 



to the classical frequency of small oscillations at the bottom of the well cu = 2SDJ1 — h 2 



while near the barrier top this frequency reduces to the instanton frequency uq. Here we 
discuss the low energy limit of the level splitting. Under the condition E <C All, the action 
W may be expanded in power series of k' according to another group of formulae 



K(jfe) = lni + - flnl - l) k' 2 + — flnl - -) k' 4 

fct A \ hi I P.A \ hi P. I 



64 
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(24) 



E W = l + iflni-i)^ + lflni-a^ 
y ' 2 V k' 2J 16 V k' 127 



Il(tanh 2 x\, k) 



1 ~ h l , l+h x , 4 l + h x + y/l-hfro^ | fe B + l fc/2ln ■ 



2/u 



-Xq + 



2/i, 



In — - 



8h 3 



8h x 



4 

Jfe 7 



'4/i? - 1 



1 — h 2 x Xo + 



3 7hl + 8hl-l\ lj4 , 9/^ + 11/^-2^4 



k'" + 



\ UKl v * u 256 h% J 128h 3 x 
and A;' 4 ~ 16E'h 2 /(l — h 2 ) 2 . Inserting them back into ( |i4|) we obtain 

2x , l-^.,,4 1-^.-- ' 



W = 2SJ1 - hi 



- 1 - 



16h 2 k' 



(25) 



'1 - h 2 x 64/>2 

which reduces to the ground state action @ at fc' = 0. This expression may be rewritten in 
a more compact form 



UOq E 



(26) 



with 



Q 



8S(l-h l) 
hi 



(27) 



Approximating the energy levels in the well by a harmonic oscillator, i.e., E — (n + |)u;o, 
and taking into account corrections from the functional-integral technique |16| , [T7| J3[1 , simplify 
(EST) into 



W = s c (0) - In 



8Se(l-hlY 



n + Z)hl 



2\ (n+i) 



(28) 



1 ) X 



so the low-lying energy level splitting takes the form flT6| , |T0 

AE, 

where the ground state splitting 



AEo 



-q' 



(29) 



AE 



8Slp 
i 

7T^ 



expx/l - hi 



2S 



h 



2S 



is proportional to power 25 of the perturbative transverse field h x WU 



(30) 



Now consider the phase transition in a biaxial anisotropic ferromagnetic grain without 
applied magnetic field 



T~t = K\S 2 Z + K 2 S% 



(31) 



which has been investigated intensively |8,13|,19|. This model possesses a XOY easy plane 
and an easy axis along the x direction with the anisotropy constants K\ > K 2 > 0. We 
show here that it provides another example which exists first-order transition from classical 
to quantum regimes. With the help of the coherent-state path integral the effective Hamil- 



tonian in terms of continuum field variable can be written in the form Ti 



2m(4>) 



and 



m(4>) 



1 



U{(j>) = K 2 S(S + 1) sin 2 



(32) 



2^(1 - A sin 2 0)' 

where A = |p < 1 and mass m{4>) is field dependent. The periodic instanton solution is 
given by jl9fl 



arcsm 



. — k 2 sn 2 (out, k) 
\ 1 — \k 2 sn 2 (cut, k) 



(33) 



with 



k 



n\-\ 



k 2 s{s + r 



E, 



(34) 



d 



where to = u y 1 — X/nf denotes the small oscillation frequency at the position of periodic 
instanton and Uq = 45(5 + l)KiK 2 . The non-monotonically decreasing behavior of the 
periodicity of the solution, t(E) = 4K(k)/u, has been pointed in Ref. [§] where the au- 
thors proved that beyond a critical value of coupling the spin system acquired a first order 
transition as a result of the field dependence of its effective mass. We now turn to evaluate 



S 



the effective free energy dependence on dimensionless energy scale P for various anisotropy 
constants. 

The Euclidean action evaluated for the periodic instanton trajectory is given by eq. ( |I~3"D 
with 



W 



K(fc) - (1 - \k 2 )Tl(\k 2 , k) . (35) 



XKt 

The second order transition temperature for this model is defined by eq.flTB]) with the in- 



stanton frequency u = y\U"(n/2)\ /m(n/2) = y/1 — Xto . Near the top of the barrier the 
third kind of elliptic integral should be expanded as 

I1(AA; 2 , k) = I + (2A + if-k 2 + (8A 2 + 4A + 3) ^-A; 4 + (l6A 3 + 8A 2 + 6A + 5) ^-ixk & . 

2 8 ^ ' 128 ^ ' 51z 

(36) 

After a straightforward calculation we obtain the expansion result for the free energy QTBj) 
in terms of k up to the sixth order 



F/AU =l-P + e(l-\)(k 2 + (6\+l)j + (40A 2 + 8A + 3) ^ . 



(37) 



The effective free energy analogous to Landau theory of phase transitions near the top of 
the barrier (P<1) reads 

F/AU = 1 + (0 - 1) P + (1 - 2A) P 2 + ^ (8A 2 - 8A + 3) P 3 + 0(P 4 ) 

(38) 

with the exact relation between k 2 and P 

k2 = T^P Y < 39 > 

The factor before P changes sign at the phase transition temperature T = T . The bound- 
ary between the first- and the second-order transition is obviously seen to be A c = |. The 
computed dependence of F on P for the entire range of energy is plotted in Fig. 2 for 
S 2 = 1000, Zfi = 1. At A = 0.3, there is only one minimum of F at the top of the barrier 
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for all T > Tq 2 \ Below it continuously shifts to the bottom as the temperature is low- 
ered. This corresponds to the second-order transition from thermal activation to thermally 
assisted tunneling. At A = 0.8, however, there can be one or two minima of F, depending 
on the temperature. The actual phase transition(in this case the first-order) occur at the 
temperature when the two minima have the same free energy, which for A = 0.8 takes place 
at T {1) = 1.122T (2) . 

The relation between tunneling at excited states and finite temperature can be under- 
stood as follows: Below the crossover temperature T the particle tunnels through the barrier 
at the most favorable energy level E(T) which goes down from the top of the barrier to the 
bottom of the potential with lowering temperature. Such a regime is called thermally as- 
sisted tunneling(TAT). The second order transition from the classical thermal activation to 
TAT is smooth, with no discontinuity of dT/dt at T , and the transition temperature is given 
by T . In the situation under which the first order transition may occur tunneling just be- 
low the top of the barrier is unfavorable, the TAT is suppressed, and the thermal activation 
competes with the ground state tunneling directly, leading to the first order transition. 

The flatness of the barrier top is not the necessary condition under which the first order 
transition may occur. For the constant mass model the more favorable potential for the first 
order transition is that the top of the barrier should be wider so that the particle doesn't 
have more advantage to tunnel through the barrier from the excited states than from the 
ground state. In our second model the pseudo-particles near the top of the barrier are 
"heavier" than those in the bottom of the well, i.e. m(ir/2) = m(0)/(l — A) > m(0). The 
advantage for the particles to tunnel through the barrier at higher excited states is again not 
very obvious, which leads to the first order transition from the thermal activation directly 
to ground state tunneling. 

A simple estimation for the crossover temperature is given by 

T<°> - AU ~ K ^ S+ ^ (40) 
Io -2S C (E = 0)- 21ri iW| ^ 

where the superscript of indicates that the ground state tunneling is considered. In Fig. 
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3 we plot the dependence of the actual crossover temperature on A. For A < |, the actual 
crossover temperature T is read out from Fig. 2(b). The escape rate can be conveniently 
represented in terms of the effective temperature defined by 

r ~ exp[--^-d = «p[-%] (41) 

The dependence T e ff(T) = AUT/F min is represented in Fig. 4 for different anisotropy 
constant ratio A. The most significant difference between the crossover temperature T and 
the actual crossover temperature T arises in the limit A — > 0, that is, /T^ = tt/A 
0. 785. 

The first-order escape-rate transition considered above is the transition from thermal 
activation to thermally assisted tunneling near the bottom of the potential and not directly 
to the ground-state tunneling || due to the logarithmic divergence of the instanton period 
t(E) for the energies near U m \ n . In some field-theoretical models, as, e.g., the reduced 
nonlinear 0(3)— a model, r approaches near the bottom of the potential. Accordingly, 
the second derivative of W(E) and F(E) is negative everywhere, as for the rectangular 
potential for particles. In such a situation, as it is clear from Fig. 2(b), the minimun 
of F(E) can only be at the barrier top or potential bottom. That is, thermal activation 



competes directly with the ground-state tunneling. It was shown that pC§ adding a small 
Skyrme term to the reduced nonlinear 0(3)— a model causes r to diverge near the bottom 
of the potential, with the accordingly small amplitude. This is, in a sense, similar to the 
situation realized in this spin model for A — > 1. 

In conclusion, we present in this letter the periodic instanton calculation for the first- 
and second-order transitions between quantum and classical regimes for two spin models. 
The level splitting formula for the excited states is derived and checked for the uniaxial spin 
system. Our results for uniaxial spin model confirm Ref. and for the biaxial anisotropic 
spin system we find that the transition from the classical regimes with lowering temperature 
is of the first order for A above the phase boundary line A c = ~, and of the second-order 
below this critical value. 
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Figure Captions: 

Fig. 1. The DSHG potential and the periodic instanton configuration for S — 10, D — 
0.6K, h x = 0.01 

Fig. 2. Effective free energy for the escape rate: (a) A = 0.3, second-order transition; 
(b) A = 0.8, first order-transition. 

Fig. 3. Dependence of the actual phase transition temperature T on the anisotropy 
constant ratio. 

Fig. 4. Dependence of the effective temperature T e ff on T for different values of A 



15 



F(P)/AU 




P=1-E/AU 



Fig. 2. (a) 




Fig. 2. (b) 
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